NOTICE:  Whea  goverEmeat  or  other  draviags>  specl- 
flcatloas  or  other  data  are  used  for  say  purpose 
other  thaa  la  coMiectloa  with  a  deflaitely  related 
goverameat  procurement  operatloa,  the  Ut 
GOverameat  thereby  lacurs  ao  reepoaslbility>  nor  any 
obligatloa  vhat soever j  and  the  fact  that  the  Goveia^ 
Sent  say  have  forsulated>  furalshed>  or  la  any  vay 
supplied  the  said  draviags,  speeifieations>  or  other 
data  is  not  to  be  regarded  by  implication  or  other** 
vise  as  la  any  manner  llcenslag  the  holder  or  any 
other  person  or  corporatioa,  or  GOnveylag  any  ri^ts 
or  penalssloa  to  manufacture/  use  or  se^  any 
patented  invention  that  may  In  any  way  be  related 
thereto. 
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ABSfRACf 

'  this  report  deals  with  the  ealculatlon  of  unsteady  lami* 
nar  houndary  layers  over  arbitrary  cylinders  in  an  incompres¬ 
sible  flow.  The  calculation  procedure  is  based  on  an  improved 
integral  solution  to  the  governing  unsteady  boundary- layer 
equations.  The  essential  feature  of  this  improved  solution 
is  to  treat  the  usual  integral  solution  as  a  first  approxjmia- 
tion,  which  is  then  used  to  linearize  the  governing  equation 
in  such  a  way  that  improved  unsteady  velocity  profiles  in  the 
boundary  layer  are  readily  obtained.  The  integral  solution 
is  first  described  jji  detaiii  based  on  two  assumed  types  of 
profiles i  One  being  the  well-loioun  fohlhausen's  fourth-degree 
polynomial 4  and  the  other >  an  exponential  function  containing 
one  arbitrary  parameter.  i,The  basis  of  the  improvement  tech¬ 
nique  is  then  presented.  Also  introduced  is  a  simple  and 
convenient  error  criterion  which  is  capable  Of  indicating# 
without  the  knowledge  of  an  exact  solution,  whether  the 
improved  solution  is  actually  more  accurate  than  the  basic 
Integral  solution.  Six  nvmierical  examples  are  then  described 
to  demonstrate  the  application  of  the  improved  integral  solu» 
tion  and  also  to  illustrate  the  use  of  the  error  criterion. 
Pinaiiyi  it  is  concluded  on  the  basis  of  the  numerical  re* 
suits  that  in  a  general  imsteady  problem  accurate  results  can 
be  obtained  by  using  the  improved  Integral  solution  based  on 
velocity  poiynorolals  within  its  whole  rsuige  of  validity,  and 
that#  outside  of  this  range#  the  integral  solution#  based  on 
exponential  fiaictlon#  should  be  used  with  or  without  the  lm« 
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provement,  depending  on  the  relative  rtagnitudes  of  the  res 
pective  error  quantities^ 


N0MENCI4TUHE 

2 

local  skin-fpictlGfi  coefflGlent,  t^/(l/tpUQ  ) 
arbitrary  function  Introduced  In  Equation  (21) 
universal  functions  defined  in  Equations  (3) 
dimensionless  stream  function  in  Equation  (39) 
form  parcaneter  in  the  exponential  profile 
function  of  t* 

ratio  of  displacement  to  momentum  thiclcnesses 
error  criterion 

vuiiversai  finiction  defined  in  Equations  (8) 
characteristic  length  of  cylinder 

function  defined  in  Equation  (29) 
fimction  defined  in  Equation  (30) 

Reynolds  number^  UpL/v 
time  variable 


velocity  components  in  the  x*  and  y-  dinectlons, 
respectively 


reference  velocity 

U«/Uq 

distance  coordinate  along  cylinder  surface  measured 
from  forward  stagnation  point  or  from  leading  edge 

x/L 

diatance  coordinate  normal  to  cylinder  surface 
measured  from  surface 


«2®/v 

ZUp/L 

unsteadiness  parameter 
V»)]  k-0 

a  boundary- layer  thicieness 
momentum  thickness 


error  quantity  defined  in  iquatlon  (3^) 


tig  y/^/W 

X  ioim  pas^aineter  In  the  polyndmial  pfdfile 

ix  dynamic  viscosity 

V  kinematic  viscosity 

i  l\inctlon  defined  in  Equations  (T) 

i*  |L/Uq 

p  density 

wall  shear  stress 
0  function  Of  X* 

Subscripts: 

l#J  grid  points  In  Fig.  l 

0  Initial  condition 

0  basit  integral  solution 

w  surface  condition 

e»  free^stream  condition 

I 


region  O^Tj^l 
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Del^ivatives  with  s^espect  to  the  Indepehdeht 
variable 


INfROplJGflON 

The  two  dlmehslOhal  Incompressible  laminar  boundary* 
layer  theory  for  steady  flows  around  cylinders  Is  now  well 
established  in  the  literature  4  it  is  unfortxmate,  however, 
that  the  corresponding  theory  for  unsteady  flows  is  much 
less  advanced,  even  thou^  it  is  generally  loio^  that  un* 
steady  flow  is  physically  more  natural  than  flow  under 
steady  conditions,  and  that  considerable  physical  insight 
to  steady“fiow  behavior  may  be  gained  by  studying  unsteady 
boundary  layers  concerning  coupled  effects  of  vorticity 
diffusion  away  from  the  cylinder  and  fluid  convection  along 
the  cylinder  surface 4  The  main  difficulty  is  evidently  the 
inclusion  of  the  time  variable  in  the  governing  laminar 
boundary^layer  equations*  Nevertheless,  many  significant 
Investigations  on  unsteady  laminar  boundary  layers  have 
appeared  in  the  literature,  especially  In  recent  years.  A 
concise  survey  on  this  subject  has  recently  been  given  by 
Stewart son  [i]  Briefly,  most  published  studies  on  vuisteady 
laminar  boundary  layers  may  be  grouped  into  six  major  areas. 
The  weli®known  classical  Rayleigh  problem  has  now  been  ex* 
tended  to  other  infinite  plate  problems  with  more  general 
boundary  conditions  [2*5] Another  area  of  study  centers 

^Numbers  in  brackets  refer  to  References  at  the  end  of  this 
report , 

^No  attempt  is  here  made  to  quote  all  References . 
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on  the  predlGtion  of  the  Instant  of  onset  of  laminar  separa** 
tion  for  unsteady  flow  starting  from  rest>  as  undertaken  by 
early  workers  In  this  field  [6]*  a  third  area  of  investiga¬ 
tion  concerns  with  exact  solutions  and  possible  exact  solu¬ 
tions  to  the  governing  differential  equations  for  a  number 
of  specific  unsteady  cases  [7*11] .  purthermore,  the  leading- 
edge  problem  of  unsteady  flow  over  a  semi- infinite  plate  has 
received  much  attention  [12-14].  Also>  very  recently  the 
effect  of  free-stream  oscillation  on  laaninar  boundary*layer 
behavior  has  been  studied  in  great  detail  by  many  Investi¬ 
gators  [15-22] .  The  last  area  of  investigation  deals  with 
general  unsteady  laminar  boundary  layers  over  arbitrary 
cylinders  [23-25] . 

In  a  receht  study  of  hydrodsnajuic  stability  of  unsteady 
laminar  boxndary  layers  over  arbitrary  cylinders,  it  has 
become  apparent  that,  similar  to  the  corresponding  steady- 
flow  problem,  the  stability  characteristic  would  depend 
strongly  on  the  second  profile  derivatives  of  the  imsteady 
velocity  profiles  in  the  boundary  layers,  and  consequently, 
accurate  determination  of  these  profile  derivatives  becomes 
necessary.  For  arbitrary  cylinders  and  arbitrary  free- 
stream  unsteadiness,  only  approximate  solutions  are  avail¬ 
able  in  the  literature,  loth  Schuh  [24]  and  Yang  [25]  have 
developed  one-parameter  integral  solutions  to  this  general 
problem  with  the  essential  difference  in  the  chosen  profiles, 
with  Schuh  [24]  utillzijig  a  combination  of  Hartree's  profiles 
for  steady  wedge  flows  and  Pohlhausen's  fourth-degree  poly- 
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nomials>  and  Yang  [25} j  profiles  derived  from  exact  sirni-* 
larity  solutions  for  \Msteady  stagnation  flows  [9].  Even 
though  these  integral  solutions  do  yield  satisfactory  over¬ 
all  unsteady  laindjiar  feoundary-layer  behaviors,  there  is 
considerable  doubt  that  these  approximate  solutions  would 
generally  give  second  profile  derivatives  accurate  enough 
for  stability  calculationSi  as  this  is  a  well  idioim  fact 
in  the  steady-flow  theory. 

This  weakness  of  the  integral  solutioni  at  least  for 
steady  cois^ressible  or  mcompressible  flows  >  can  be  corrected 
to  a  very  significant  extent  by  an  toprovement  technique 
developed  recently  [t6,27l .  since  the  vinderlying  idea  of 
this  technique  applies  equally  well  to  the  general  unsteady 
lafflinar  boundary^layer  problem,  there  is  good  indication 
that  it  may  be  directly  extended  to  the  present  problem  to 
obtain  improved  unsteady  velocity  profiles.  However,  in 
view  of  the  complexity  introduced  by  the  added  time  variable, 
such  expectation  can  only  be  ascertained  by  detailed  analysis 
of  the  unsteady  problem  and  con^arison  of  calculated  results 
with  loiown  accurate  solutions  in  the  literature.  The  primary 
pu^^ose  of  this  report  is  to  present  such  an  analysis  and  the 
corresponding  results,  and  in  addition,  it  will  also  be  shown 
that  the  general  validity  of  this  improved  integral  solution 
may  be  readily  determined  by  an  error  criterion  [27]  . 
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PORMULAflON  Ai©  INTEGRAL  SOLUTION 

When  an  arbitrary  two-diinensional  cylinder  moves  \insteadi* 
ly  in  an  inGompressible  stagnant  fluid,  unsteady  laminar  boundary 
layer  develops  on  the  cylinder  surface,  in  the  incompressible- 
flow  theory,  this  problem  is  exactly  equivalent  to  one  of  the 
same  fluid  flowing  over  the  stationary  cylinder  with  the  same 
unsteady  velocity.  Thus,  with  the  coordinate  system  fixed  on 
the  cylinderj  the  well-ioioiim  laminar  boundary-layer  equations 
may  be  written  as  follows: 

du  ^  ^  d-u 

+  u^  +  v-^  =  •  +  u«?— ^  ( 1 ) 

dt  dx  dy  dt  dx  dy^ 


du  dv 

^  =  0  (2) 
dx  dy 

where  x  is  the  coordinate  along  the  cylinder  surface  measured 
from  the  forward  stagnation  point  or  the  leading  edge,  y 
the  coordinate  normal  to  the  cylinder  measured  from  the  sur¬ 
face,  u  and  V,  velocity  components  in  the  x-  and  y- 
directions,  respectively,  t  the  time  variable,  v  the  kinematic 
viscosity,  and  the  subscript  w  indicates  conditions  in  the 
free  stream.  The  Initial  condition  generally  depends  on  the 
Specific  problem,  while  the  boinidary  conditions  may  be  written 


y-0  U^v=0 

u— ^^u^(x,t) 


(3) 


y— 


For  arbitrary,  but  pfescrlbed  variations  of  u^(x,t)>  Integral 
solutions  based  on  somewhat  different  assumed  velocity  pro^ 
files  have  been  presented  by  schuh  lJ4]:  and  Yang  [25]  •  The 
present  analysis  follows  essentially  the  latter >  except  for 
some  details,  as  will  be  shovm  later  in  the  report.  When  the 
momentum  equation  (1),  after  eliminating  v  according  to  (2)> 
is  integrated  across  the  boundary  layer,  the  following  well> 
knowi  integral  equation  results: 


■t  (2  H) 


V  dx 


\.^2 


W 


where  is  the  usual  momentum  thickness >  H  the  ratio  of 
displacement  to  momentum  thicknesses,  the  wall  shear 
stress,  and  u  the  viscosity*  Now  if  we  assume  that  the  un¬ 
steady  velocity  can  be  represented  by  a  one»parameter  femily 
of  curves,  i.e. 

u*  =  u*(x,Ti)  (3) 

^00 

where  t]  =  y/6,  x  *  x(6,u^),  and  6  is  an  unsteady  boundary-^ 
layer  thicloness,  Equation  (4)  becomes  a  first-order  partial 
differential  equation  for  the  unioiown  6(x,t).  In  order 
to  facilitate  the  solution  to  this  equation.  It  may  be  con¬ 
veniently  cast  in  a  simpler  dimensionless  form  as  follows  {2^ 
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^  dK 

U«  +  F-i 

dx*  ^ 

H 

/O  I  /O 

*‘h 

K(C 

«  A 

+  Hiy 

(6) 

Where 

K  »  zi 

Z  « 

V 

i  ^ 

^00 

hx 

}  ^00 

u  dt 

00 

R 

H 

+.*  „  „o_ 

L 

111 

’  % 

u«  * 

^0 

(7) 

4  35* 

A  rt 

4*  dx* 

B  « 

1  dt* 
dp 

c  - 

^  dp 

and  Uq  Is  a  refefsnce  velocity,  and  L>  a  characteristle 
length  of  the  cylinder*  Rirtheinnore,  universal  functions 
Kj  Fg  and  F^  are  pure  ftmctions  of  the  form  parameter 
x,  which  may  be  written  as 


K  > 


M-U«  d 


dM  dX 

F,  •  H 

*  dX  dK 


Fg  *  2F2  ^  aiffl 


(8) 


and  In  addition,  H  and  (dg/d)  are  also  wilversal  functions 
of  the  same  parimeter  X,  provided  that  x  is  sptolflcally 


defined  as 


V 


(9) 


which  obviously  reduces  to  the  Pohlhausen's  form  parameter 
for  steady  flows*  It  now  becomes  clear  that  even  though  K 
Is  explicitly  written  as  the  unlmown  In  iquatlen  (€),  the 
form  parameter  x  is  actually  the  unloiown  to  be  sought j  In 


view  of  the  unique  delation  (8)  between  K  and  X*  once  the 
velocity  profile  (5)  is  asstmed. 

iefore  solutions  to  Equation  (6)  can  be  attempted >  it  is 
necessary  first  to  introduce  specific  profiles  in  the  form 
of  Equation  (5)  ^d  to  determine  the  corresponding  universal 
functions.  As  already  pointed  out  previously,  Schuh's 
analysis  £24^  utilizes  both  the  Pohlhausen's  fourth^degree 
velocity  polynomial  (for  K40)  and  Hartree's  profiles  for 
steady  wedge  flows  (for  K>o),  and  Yang  [253  based  his 
choice  on  the  exact  profiles  for  unsteady  stagnation  "hyper¬ 
bolic  time*variation"  flows  [9].  I'he  obvious  wealcness  in 
schuh's  choice  Is  that  the  resulting  universal  functions  are 
discontinuous  at  K  *  0,  and  this  point  of  discontinuity  is 
expected  to  occur  quite  frequently  in  any  general  unsteady- 
flow  problem.  In  addition,  the  Hartree's  wedge-fiow  profiles, 
as  well  as  the  corresponding  universal  functions,  are  only 
known  In  tabulated  forms,  and  hence  cannot  be  manipulated 
readily.  Yang's  profiles,  though  continuous  throughout  the 
entire  region  of  K,  evidently  suffer  the  same  difficulty. 
However,  since  the  present  improved  Integral  solution,  as 
will  be  described  in  the  next  section,  utilizes  the  integral 
solution  as  a  first  approximation,  the  velocity  profiles  and 
the  associated  universal  functions  from  the  integral  solution 
must  be  manipulated.  Consequently,  all  these  available  un¬ 
steady  profiles  are  not  suitable  for  the  present  use.  In 
the  present  study,  two  distinctly  different  profiles  have 
been  considered,  One  is  still  the  Pohlhausen's  fourth-degree 
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pdlyndmlal>  whichi  hdwevei*>  Is  alldwed  td  be  valid  In  the 
range  12,  with  the  Idwer  limit  Gdrrespdnding  te 

laminar  separatldn.  Fdr  x>12,  the  veldclty  pdlyndmlal  be* 
cdmes  ndn*senslblei  and  hence  is  td  be  replaced  by  an  ex* 
pdnential  prdfile.  This  discdntinuity  is  net  expected  td  be 
tdd  serldus  in  actual  applicatidns,  since  it  dnly  dccurs  in 
a  regidn  df  large  values  df  cdrrespdnding  td  very  large 
acceleratidn  in  the  free  stream*  Ndw  the  details  df  these 


twd  assumed  prcfiles  are  separately  described. 

The  Pdhlhausen's  fdurth  degree  pdlyndmlal  fdr  the  velccity 
prdfile  Is  well  kndwn,  and  may  be  written  as 


frdm  which  imiversal  functidns  in 
mined j  resulting  in 


o^Ti^l 


t)  >  1 

($)  can  be  readily  deter* 


3  X 
10  120 


BT 

X 

313 

9^5  9072 

X 

37  X 

(2  +  ^  *  — ) 
6  315  9^5  9072 


X^  ^  X  i3X^  X^ 

9072  1375  ilB^OO  272160 


(11) 
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Outside  the  region  of  validity  of  the  above  poljhiomial. 
Equation  (10)  Is  to  be  replaced  by  an  exponential  function, 
which  has  been  utilized  in  both  steady^  and  unsteady^flow 
cases  [28,29,3]  with  high  degree  of  success  especially  in 
the  unsteady^flow  problems^  It  is  in  the  following  form: 

u*  =  1  =  -  GT1)  0471^«b  (12) 


where  ti  =  y/d  and  the  form  parameter  0  plays  the  same 


role  as  x  in  the  pohlhausen's  polynomials  It  should  be 
emphasized  here  that  the  boundary»layer  thickness  6  In 
Equation  (It)  has  no  relation  whatever  to  that  In  Equation 
(10).  In  addition,  physically  sensible  profiles  are  only 
obtainable  for  negative  0  values,  with  a  =  -^l  for  laminar 
separation.  It  is  also  seen  that  boundary  conditions  at  ln» 
finity  are  automatically  satisfied,  and  the  relationship  be^ 
tween  0  and  5  can  be  easily  obtained  from  evaluating  the 
momentum  equation  (l)  at  the  wall,  l.e. 

1  +  20  ^  —4  (13) 

V 


Now  once  the  profile  Is  knoum,  the  corresponding  universal 
functions  may  again  be  determined  explicitly.  They  arc  shown 


as  follows: 

K  » 


(1  »»■  2G)( 


1 

2 


G-  2  ^2  '  1  ^ 

2  6  2  2  4 


^(1  -  Q) 

2  -  2G  -  G^ 


Pg  ^  j(l  +  G)(2  i.  2G  -  0-) 


(14) 


P, 


12 

8  +  50 


/ 


I 


0^(2  -  2G  -  G^) 
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Now  solution  to  the  unsteady  integral  equation  is  oon- 
sidered*  Equation  (6)  is  a  linear  first-order  partial  dif¬ 
ferential  equationi  the  most  general  solution  of  whiGh  is  by 
the  method  of  eharaGteristies .  However >  this  method  yields 
exact  sdiutions  to  Equation  (6)  only  in  those  cases  where 
two  integrable  equations  can  be  obtained  from  the  correspond¬ 
ing  characteristic  equations*  For  general  unsteady  free- 
stream  velocity  distribution,  other  means  of  solution  must 
be  used  instead,  in  view  of  the  above,  problems  of  special 
interest  may  be  grouped  Into  the  following  classes,  relative 
to  different  geometry  and  methods  of  solution: 

1 .  stagnation  Plow, with  Free^gt ream  Velocity  Varying  Arbitrarily 
with  Time 

For  this  class  of  problems,  the  free-stream  velooity 
distribution  may  be  written  as 

where  g  is  an  arbitrary,  but  continuously  differentiable 
function  of  time.  With  I*,  A  and  B  given  by  pure  functions 
of  t  and  C  =  O,  Equation  (6)  reduces  to 

dJC  I*  r 

—  +  K(ffiB-A)J  (15) 

dt  J'i  '  ' 

Which  may  be  conveniently  solved  nufflerlcally  for  the  untaiown 
X  or  0,  depending  on  the  profiles  chosen.  In  terms  of  cer-* 
tain  initial  condition.  For  instance,  if  g(0)  =?=  I,  corres* 
ponding  to  initially  steady  motion,  then  X(Q)  and  0(Q)  can 
be  determined  from  ate®<iy*'Flw  theory. 


2.  Flow  Qver^an_Inflnlte,  Plate, jiflfch_Pree?iSj£p^^ 

Varying  Arbitrarily  with  Time 

With  =  g(t*)  and  all  derivatives  with  respect  to 
X*  neglected >  It  may  be  readily  shown  that  iquation  (6) 


now  as  Slimes  the  following  form: 

,* 


dK 

IP 


£ 

F. 


F^  +  Biflc] 


Where  both  4*  and  B  are  now  piije  functions  of  t**  The 
above  equation  may  again  be  numerically  integrated  in  terms 
of  any  specific  initial  condition  depending  on  the  physical 
problem  under  consideratlon>  resulting  in  either  x(t*)  or 
Q(t*)*  It  may  be  pertinent  to  mention  here  that  the  corres* 
ponding  problem  of  unsteady  flow  over  a  semi* infinite  plate 
presents  some  difficulty  regarding  its  solution,  it  must  be 
solved  in  the  same  marmer  as  a  general  problem  of  an  arbitrary 
cylinder  in  an  arbitrarily  imsteady  flow,  as  will  be  described 
later.  One  exception  is  the  case  of  a  seml*inflnlte  plate 
moving  from  rest  with  a  constant  acceleration.  The  corres* 
ponding  integral  equation  (6)  has  a  siroilerity  solution 
satisfying  the  following  ordinary  equation  [24,25] 


dK  2P2  *  3KH 
df  1  «  2PP^ 


where  ?  ^  xV(lopt*),  which  is  subjected  to  the  boundary  con» 
ditlons  K(0)  s  0  and  K(«)  ^  K^,  where  the  constant 
is  determined  from 


^  19  * 


Since  only  one  feoundary  condition  is  necessary>  solution  to 
Equation  (l?)  is  obtained  by  matching  the  two  solutions  sat“ 
isfying  the  two  boundary  conditions  at  a  proper  value  of  p. 
other  details  of  this  problem  will  be  shown  in  a  later  sec* 
tion.  Another  exception  is  for  flow  over  a  semi*infinite 
plate  with  step*change  in  free»stream  velocity*  since  for 
this  problem  K  is  identically  zero>  all  universal  functions 


become  constants 4  and  Equation  (6)  In  this  limit  reduces  to 
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2Fc 


where  z  =  Zu^/L*  This  equation  may  be  conveniently  solved 
by  the  method  of  characteristics.  The  corresponding  char* 
acteristic  equations  are  then  given  by 


dx*  dt*  dZ* 


(20) 


from  which  two  Integrable  equations  may  be  readily  obtained. 
After  eliminating  the  two  arbitrary  constants  of  integration^ 
we  obtain  the  general  solution 
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where  the  arbitrary  function  P  is  deteralned  from  the  Inl* 
tlal  and  boundary  conditions.  For  a  flow  starting  from  rest* 
the  conditions  Z*(x*,0)  =  O  and  Z*(0,t*)  ^  0  must  all  be 
satisfied.  Hence  Equation  (21)  becomes 
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where  l^uatlofi  (22a)  indicates  the  transient  behavior >  while 
Equation  (22b )i  the  steady^state  solution,  when  the  flow 
changes  from  an  originally  steady^state  condition,  we  have 
accordingly  z*(x*,o)  *  Z^*(x*)  and  z*(o,t*)  =0.  A  simi- 
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where  Uojji^  is  the  initial  steady  free-stream  velocity#  which 
may  be  taken  as  unity  without  loss  of  generality. 

3.  Arbitrarily  unsteady  Flow  over  Arbitrary  Cylinders 

■The  validity  of  the  method  of  characteristics  as  applied 
to  this  general  problem  has  been  implicitly  assumed  by  $chuh 


[2^1 .  However#  limitations  to  this  method  were  not  pointed 
out,  nor  was  any  example  illustrating  this  application  to  a 
general  unsteady  problem  given.  Yang  [2$]  later  has  proposed 
a  step"by»step  calsnlation  procedure#  also  based  on  the  method 
of  characteristics#  to  this  general  problem.  Unfortunately# 
a  close  examination  of  this  procedure  reveals  that  it  yields 
accurate  results  only  in  cases  where  the  quantity  dK/dx* 
is  a  weak  function  of  t  .  In  the  present  study,  a  new  cal^ 
culation  procedure  is  introduced,  which  has  the  capability  of 
obtaining  numerical  solution  to  Equation  (6)  as  accurately 


as  desired.  This  procedure  is  based  on  a  lumped  approximation 
of  tho  vsrintion  of  K  with  respect  to  x  •  If  we  now  con^ 


sidef  the  grid  system  in  the  x*  *  t*  plane  as  sho^  In  Pig, 
li  It  is  seen  that  the  derivative  ^K/^x*  at  any  point  may 
fee  conveniently  approximated  fey 


/„)  ^  ^i»i>i  "  %-ii  j 

dx*  *  „  2AX* 

t  ^  jht* 

Which  involves  an  error  of  the  order  of  (Ax  ) -  , 
equation  (6)  may  now  fee  written  as 


(24) 


The  integral 


dK  i*  r  n  n*  dK 

[fj.KIO.A,®)]  .-(p)  (25) 

Por  a  given  profelem,  some  initial  condition  K(x*>o)  must 
fee  known*  Por  a  flow  starting  from  rest>  K(x*>o)  m  k^, 
which  must  %e  deter®ined  from  a  limiting  process  similar  to 
that  of  ofetaining  Equation  (18)  *  When  the  flow  is  ini*- 

tialiy  steadyi  k(x*#o)  may  then  fee  evaluated  from  the  steady* 
flow  theory,  iefore  solution  to  gquatlon  (25)  is  considered, 
it  IS  first  necessary  to  ofetain  K(Q,t*)*  For  a  felunt*nosed 
cylinder,  K(0,t*)  is  ofetalned  fey  directly  integrating  Equa¬ 
tion  (25)  with  Uqb  -  0,  which  Is  evidently  Identical  to  Equa* 
tlon  (15).  When  the  cylinder  has  a  sharp  leading  edge, 

K(o,t*)  feecomes  identically  zero.  Once  K(x*,o)  and 
K(0,t*)  are  known,  solution  to  Equation  (25)  may  then  pro* 
ceed  as  follows^?  Since  at  t*  =  Q  the  derivative  dK/dx* 
is  known  exactly,  there  Is  no  need  to  evaluate  this  quantity 
from  Equation  (24).  Prom  the  given  unsteady  free-stream 


9 

-Xn  actual  applications,  it  is  deslrafeie  to  treat  either  x 
or  Of  as  the  dependent  variafele,  instead  ef  K,  in  view  of 
the  fact  that  all  imlversal  functions  are  explicit  functions 
of  X  or  0. 


velocity  disti^^ilsution,  quantities  such  as  4  >  A>  C  and 
u^  may  then  be  evaluated  at  this  instant.  Consequently > 
now  the  night*hand  side  of  Equation  (23)  can  be  determined 
accordingly  for  all  points  along  the  x*“a3cis  (Pig.  1)* 

Then  one  step  of  numerical  integration  with  respect  to  t* 
at  each  of  these  points  immediately  yields  K(x*>nt*).  in 
this  regard i  it  has  been  fotuid  that  the  usual  Runge^Kutta 
integration  routlnei  which  is  accurate  to  the  order  of  (nt*)^ 
is  very  satisfactoryi  Together  with  the  imown  value  of 
K(OiAt*)>  derivatives  (Sk/^x*)  at  t*  =  At*  may  now  be 
evaluated  for  all  points  along  the  x*“axis  according  to 
(24).  ^is  process  is  then  repeated  for  successive  time  in* 
stints  until  a  desired  time  elapse  is  covered,  it  is  noted 
that  this  step*by-etep  calculation  procedure  has  one  advantage 
over  that  of  pure  finite*difference  numerical  approximations 
in  both  X*  and  t*,  and  that  is  that  here  the  question  of 
stability  of  the  numerical  solution  never  arises.  Purthermore 
the  present  procedure  can  be  conceivably  improved  in  Its 
accuracy  by  using  a  higher*order  differentiation  formula  than 
that  in  Equation  (24). 

Once  the  parameter  K  becomes  a  known  function  of  x* 
and  t*>  the  corresponding  fonn  parameter  X  or  0  can  be 
evaluated#  yielding  6  =  6(x*,t*).  Equation  (10)  or  (12) 
then  immediately  gives  the  unsteady  velocity  profiles.  Thus, 
the  integral  solution  may  now  be  considered  as  complete. 


AN  IMPROVED  INfEGRAL  SOLUTION  AND  A  RELATED  ERROR  CRITERION 
As  generally  recognized  In  the  application  Of  Karman* 
Pohlhausen's  integral  procedure,  there  are  two  essential 
weaknesses  which  give  rise  to  the  Inaccuracy  of  the  approxi* 
mate  solution*  One  is  that  the  governing  partial  different 
tial  equation  is  only  satisfied  -m  the  mean,  and  the  other, 
that  the  choice  of  assumed  profiles  is  entirely  arhltrary, 
and  the  resulting  accuracies  of  different  types  of  profiles 
could  be  different  by  significant  amounts*  in  the  laminar 
boundary* layer  theory,  most  earlier  improvements  for  general 
problems  have  been  based  either  on  the  choice  of  novel  pro* 
files  or  on  the  use  of  additional  integral  equations,  or  some* 
times  known  as  moment  equations.  Unfortunately,  none  of  these 
improvements  has  general  validity  to  arbitrary  problems.  This 
has  led  the  present  writer  to  develop  an  essentially  different 
Improvement  technique  for  the  general  laminar  boundary^layer 
problem  [26]  in  an  attempt  to  correct  significantly  the 
aforementioned  weaknesses  of  the  basic  integral  solution. 

This  technique  basically  utilizes  the  integral  solution  as  a 
first  approximation,  which,  when  substituted  into  the  govern* 
ing  differential  equation,  linearizes  the  equation,  which  may 
then  be  readily  solved  for  the  improved  profiles  in  closed 
forms.  In  view  of  its  apparent  success,  this  technique  has 
since  been  extended  to  other  physical  problems  [30,31]  . 

Very  recently,  it  has  been  further  improved  in  a  minor  detail 
in  that  now  the  boundary  condition  at  infinity  is  exactly 
satisfied,  and  in  addition  a  simple  error  criterion  is  in* 
troduced  such  that  validity  of  the  improvement  can  be  readily 


determined  without  Icriowledge  of  the  exact  solution  [^7]  . 

As  mentioned  previously,  this  improvement  technique  has  only 
SO  far  heen  developed  for  a  paraboiic^type  of  partial  differ* 
ential  equations  in  two  independent  variables,  its  extension 
to  equations  with  three  independent  variables  raises  some 
uncertainties i  fhe  primary  purpose  of  this  section  is  to 
present  i;ne  detailed  analysis  of  this  extension  to  the  general 
unsteady  boundary*layer  problem  in  two*dimensional  Incompres* 
sible  flow.  To  Simplify  the  presentation,  the  following  will 
be  described  only  for  the  integral  solution  based  on  pohl* 
hausen's  velocity  polynomial.  However,  the  formulation  of 
the  improved  procedure  may  be  easily  extended  to  cases  where 
exponential  profiles  (12)  are  used,  in  fact,  it  is  only 
necessary  to  replace  the  form  parameter  x  by  the  quantity 
(1+20)  in  view  of  Equation  (13)  and  the  usual  definition 
of  X.  Neverthless,  significant  differences  will  be  pointed 
out,  whenever  desirable. 

Prom  the  basic  Integral  solution  as  described  previously, 
both  X  and  6  become  known  functions  of  x*  and  t*.  in 
order  to  facilitate  the  linearization  of  the  goyeming  differ* 
ential  equations  by  introducing  the  integral  solution.  It  is 
now  desirable  to  transform  these  equations  from  the  (x,y,t) 
coordinate  system  to  the  (X,T),t  )  system.  This  transforma* 
tion  is  govexned  by 
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The  continuity  equation  (2)  Is  first  considered.  Solving  for  v 
and  transforming  i  we  have 
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(27) 


Now  the  momentwft  equation  (1)  may  be  similarly  transformed. 
With  the  help  of  Equation  (27)  and  the  following  identities: 
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Equation  (l)  leads 4  after  some  m^ipulations^  to 
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It  Is  rioted  that  this  tr^sformation  Is  exact  i  arid  the 
oplglnal  rioh^llnearitles  ape  riow  iiicluded  in  arid  Pg. 
Cdrisequentlyi  any  solution  that  satisfies  Equation  (28) 
is  also  a  solution  to  the  tnomentuin  equation  (l)^  in  the 
ppesent  linppoveraent  technique.  Equation  (28)  is  fipst 
lineapized  hy  evaluating  p^  and  Pg  on  the  hasis  of  the 
basic  integpal  solution,  and  then  integPated  to  give  pe^ 
fined  profiles.  The  Justification  comes  from  the  fact  that 
P|  and  Pg  in  Equations  (29)  and  (30)  only  involve  pro»> 
files  and  their  integrals,  and  hence  the  basic  integral 
solution  is  expected  to  describe  these  functions  rather 
accurately.  The  inaccuracy  of  the  integral  solution  in 
predicting  profile  derivatives,  especially  the  second  de* 
rivatlves.  Is  now  significantly  reduced  by  satisfying  the 
governing  differential  equation  (1)  or  (28)  much  more  close* 
ly.  Therefore,  Equation  (28)  may  now  be  approximately 
written  as 


where  P^q  attd  pg^  are  P^  and  Pg  evaluated  froffl  the 
basic  integral  solution,  respectively.  The  corresponding 
boundary  conditions  are,  for  any  x  and  t* 

T1  =  0  u*  »  0  *  u*“^  1  (32) 


In  general,  solution  to  the  linear  equation  (3I)  c^  be 
expressed  in  closed  form.  However,  in  view  of  the  composite 
nature  of  the  assumed  profile  (10)  in  the  integral  solution, 
it  is  necessary  to  integrate  Equation  (3I)  separately  in 
the  two  regions  O^TI^l  and  resulting  in 
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Where  the  subscript  i  and  II  refer  to  the  regions  0$T|5i 
and  14T1^«,  respectively,  and  Cj^,  Cg,  C-^  and  are 
constants  of  Integration^  syad  generally  functions  of  x  and 
t*.  In  order  to  eliminate  these  constants,  two  conditions 


are  required  In  addition  to  those  in  Equations  (32),  and 
they  are  the  matching  conditions  at  T]  =  1,  given  by 
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Since  both  P^^q  and  P20  ^?e  continuous  throughout  the 
region  of  T)  ssid  Equations  (34)  insure  the  continuity  of 


the  profile  and  Its  first  derivative,  it  is  seen  that  the 
new  profile  (33)  is  at  least  continuous  in  Its  second  de* 
rivatlve  in  view  of  equation  (31).  It  is  pertinent  here 
to  mention  that,  when  the  exponential  profile  (13)  is 
assumed  in  the  haslc  integral  solution,  the  question  of 
this  matehlng  in  profiles  does  not  arise >  and  it  is  only 
necessary  to  extend  the  integrals  in  (33)  to  all  values  of 
Tj  and  evaluate  the  two  constants  of  integration  aecording 
to  conditions  in  (33). 

As  mentioned  previously,  many  improvements  of  the  haslc 
integral  procedure  have  appeared  in  the  literature  over  the 
years*  Almost  without  exception>  these  improvements  are 
Justified  by  comparing  results  of  both  the  improved  solutions 
and  the  hasic  integral  solutions  with  icnown  exact  solutions 
in  certain  specific  cases.  One  major  shortcoming  of  such 
Justifications  is  that  the  accuracy  of  these  improvements 
is  never  too  well  ascertained  in  cases  for  which  no  exact 
solutions  are  available.  At  least  for  steady'S^flow  theory, 
this  shortcoming  has  now  been  corrected  by  using  a  simple 
error  criterion,  as  shown  very  recently  by  the  present  writer 
^27]*  gince  there  still  exists  some  uncertainty  as  to  whether 
or  not  the  present  formulation  of  the  improvement  technique 
does  improve  the  result  based  on  the  integral  solution 
alone  for  the  general  unsteudy'-flQw  problem,  the  use  of  such 
an  error  criterion  is  indeed  extremely  desirable.  It  is 
now  to  be  shown  that  the  development  of  such  an  error  criterion 
for  the  present  problem  presents  no  added  complexity  as 


compared  to  that  of  the  steady* flow  theory* 

Since  the  transformation  from  Equation  (1)  to  Equation 
(28)  is  exact,  it  is  clear  that  the  accuracies  of  the  in* 
tegral  solution  with  and  without  the  present  improvement 
depend  largely  on  how  closely  Equation  (28)  is  approximated 
in  eaoh  case*  Consequently,  we  may  define  an  error  quantity 
€  by 

„  du* 

)  .  -p.  ♦  Pj  -s-  -  Pg  (39) 

which  is  identically  zero  for  an  exact  solution*  For  a  given 
assumed  profile  in  the  basic  integral  solution,  deviations 
from  zero  for  e  Indicate  levels  of  inaccuracy  of  the 
approximate  solution*  in  view  of  the  fact  that  the  error 
quantity  e  defined  in  Equation  (35)  could  have  both  positive 
and  negative  values  throughout  the  entire  region  of  Ti> 
it  is  more  convenient  to  consider  a  directly  related  error 
quantity  defined  as  follows: 
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where  H  is  characteristic  Reynolds  number p  u^L/v,  It  is 
noted  that  this  definition  of  J  differs  only  slightly  from 
that  of  the  steady^fiow  theory  [27]  in  the  Inclusion  of  the 
boundary^iayer  thickness  6,  and  represents  an  average  error 
of  the  profile  details.  To  determine  the  relative  degrees 
of  accuracy  of  the  integral  solutions  with  and  without  the 
present  improvement j  it  is  then  only  necessary  to  determine 
and  compare  respective  3  values  at  different  combinations 
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of  X  and  t  i  Evidentlyj  smaller  values  of  J  correspond 
to  the  more  accurate  profiles  *  in  the  actual  determination  of 
€,  Equation  (35)  may  be  used  directly  for  the  basic  integral 
solution*  For  the  improvement,  a  further  slmpllfiGatlon  is 
possible*  since  the  improved  profile  satisfies  Equation  (31), 

It  may  be  readily  shown  that 

au* 
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Finally,  it  is  particularly  important  to  observe  that  compari* 
sons  can  only  be  made  of  J  values  based  on  the  same  type  of 
assiimed  profiles  used  in  the  basic  integral  solution*  Differ* 
enl  types  of  initially  assumed  profiles,  such  as  polynomials 
and  exponential  functions,  inevitably  lead  to  d's  which  have 
entirely  different  meanings,  since  these  6  quantities  are 
directly  involved  in  the  transformation  which  eventually  leads 
to  Equation  (35)j  meaningful  comparisons  of  the  corresponding 
J  values  evidently  caamot  be  expected » 

The  Improved  Integral  solution  presented  in  this  section 
specifically  deals  with  the  improvement  of  the  local  unsteady 
velocity  profiles.  In  particular,  the  boundary»"layer  thlcigiess 
function  6  remains  the  same  as  that  given  by  the  basic  Integral 
solution.  Obviously,  this  is  done  for  the  simple  reason  to 
keep  the  calculations  at  a  reasonable  level*  Nevertheless,  it 
is  also  conceivable  that  a  similar  foimulation  as  that  used  in 
the  steady*»flow  theory  ^27]  can  be  carried  out  to  successively 
approximate  solution  to  Equation  (28)  with  an  iteration  scheme 
based  on  the  integral  solution  as  the  zeroth*order  approxlma* 


tion>  with  the  view  of  possibly  obtaining  an  exact  solution  * 
However >  in  view  of  the  greatly  increased  complexity  of  the 
unsteadyBflow  problem,  such  ^  iteration  solution  does  not 
seem  to  be  justifiable. 

NOI^ICAL  EXAMPLES 

For  the  purpose  of  demonstrating  the  application  of  the 
improved  integral  solution,  as  well  as  illustrating  the  use  of 
the  error  criterion,  several  numerical  examples  have  been  cal* 
culated  in  detail  on  an  iSM  16^0  digital  computer.  Wieaever 
possible,  both  assumed  profiles,  polynomial  and  exponential 
function,  are  used  in  the  basic  integral  solutions  such  that 
range  of  validity  of  each  of  these  profiles  may  be  determined . 
These  examples  and  their  results  are  now  described  individually. 
EXAMPy:  i.  Stagnation,  Hyperbolic  Time? Variation  Flow 

Similarity  solutions  to  the  unsteady  two*dimensional  in* 
compressible  boundary*layer  equations  have  been  considered  by 
many  investigators  [J*ll]  .  However,  only  a  few  detailed  exact 
velocity  profiles  are  loiowi.  Yang  [9}  has  considered  the 
free*streara  velocity  variation  given  by 


where  d  is  an  unsteadiness  parameter,  positive  for  accelera* 
tion  and  negative  for  deceleration.  This  is  loiown  as  the 
stagnation  hyperbolic  time*variation  flow  for  which  a  simi* 
larity  solution  exists.  The  associated  ordinary  differential 
equation  assumes  the  following  form: 
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-  (f‘)^  +  i  =  a(»*  f"  +  f‘  «  1)  (39) 

a 

which  1$  subjected  to  the  boundary  conditions 

f(0)  =  f'(0)  «  0  f '(~)  ^  1 

where  f '  is  the  velocity  profile  u*  and  prime  denotes  de*^ 
rivatives  with  respect  to  the  similarity  variable  ^  defined 
as  y(ueo/vx)  i  Numerical  solutions  to  Equation  (39)  have 
also  bean  given  by  yang  [9]  for  a  wide  range  of  a  from  *3*o 
to  +l*6i  The  available  exact  velocity  profiles >  and  in  parti* 
cular  the  profile  derivatives  provide  us  with  an  excellent 
opportunity  to  assess  the  validity  of  the  present  improved  in* 
tegral  solutions  based  on  both  types  of  assumed  profiles. 

For  this  problem>  it  may  be  readily  shown  that  A  =  Vd+n)* 
B  =  a/(l+a)  and  C  »  0.  Since  K  is  now  a  constant,  dependent 
only  on  o,  dK/dt*  becomes  identicaily  zero.  Hence  Equation 
(15)  reduces  to 

(1  +  a)p^  +  K(aH  *  A)  »  0  (40) 

which  immediately  yields  x(a)  and  G(a).  When  these  are  sub* 
stituted  in  Equations  (10)  and  (12),  respectively,  velocity 
profiles  Uq*  are  obtained.  This  then  completes  the  integral 
solution.  For  the  improvement  solution,  functions  and 
Pg  in  Equations  (29)  and  (30)  now  become 
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For  the  exponential  function^  it  is  only  necessary  to  replace 
X  by  (1  +  2G)»  These  equations  are  not  only  used  to  evaluate 
P|Q  and  PgQ  from  the  integral  solution^  but  also  to  deter¬ 
mine  P|  and  Pg  from  the  improved  profiles  such  that  the 
error  quantity  e,  and  eventually  J  may  be  calculated  according 
to  Equations  (37)  and  (36) ^  respectively.  Based  on  both  fourth- 
degree  polynomial  and  the  exponential  function^  the  integral 
solutions  with  and  without  the  present  improvement  have  been 
obtained  by  detailed  calculations  for  a  values  of  1,6,  l*o> 

0>  -1.0,  -2,4  and  -2.8*  All  corresponding  J  values  have 
also  been  evaluated,  according  to 


2(1  +  a)  1  +  a 


u  an 


u*(u*  ^  a) 


X  4.  ry. 
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(polynomial) 


(exponential) 


fable  1  shows  some  results  of  these  calculations  in  terms  of 
the  surface  derivatives  of  the  velocity  profiles  f"(0) 
and  the  error  quantities  J//i  »  at*,  integral  solutions  based 
on  polynomials  only  exist  up  to  an  a- value  of  0,2,  while  those 
based  on  the  exponential  function  does  not  seem  to  have  an 
upper  limit.  Lower  limits  for  both  assumed  profiles  obviously 
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cofresporid  to  laminar  separation.  The  results  shown  in 
fahle  1  are  very  instructive  In  several  respects  *  Even 
though  the  error  quantities  only  indicate  an  average  for  the 
entire  profiles  and  yet  f''(o)  values  refer  to  only  ohe 
specific  value  of  Ti,  there  is  still  a  good  correlation  be* 
tween  the  two^  namely >  lower  values  of  the  error  quantity 
do  correspond  to  closer  agreements  with  the  exact  solution ^ 
Apparently,  these  surface  profile  derivatives  are  represen* 
tative  of  other  local  profile  values  as  far  as  accuracy  is 
concerned.  For  the  case  of  assuffied  velocity  polynomial, 
it  is  seen  that  the  present  improvement  reduces  the  error 
in  the  basic  integral  solution  considerably,  indicating  a 
corresponding  increase  in  accuracy  of  the  result*  This  de* 
finitely  suggests  the  validity  of  the  improved  solution  for 
general  use  within  the  entire  range  of  validity  of  the  basic 
Integral  solution  based  on  foh] hausen's  velocity  polynomials. 
However,  such  a  general  conclusion  cannot  be  made  for  the  in* 
tegral  solution  based  on  exponential  profiles.  For  negative 
Values  of  a,  the  present  solution  indeed  represents  improve* 
ment.  However,  when  u  becomes  positive >  J*vaiues  for  the 
present  solution  actually  exceed  that  of  the  basic  integral 
solution,  indicating  that  no  Improvements  are  realized  in 
these  cases,  consequently,  the  present  solution  should  not 
be  used,  it  is  emphasized  here  again  that  devalues  for  the 
two  assumed  profiles,  which  are  seen  to  be  different  by  some 
orders  of  magnitude,  bear  no  relation  to  one  another,  due  to 
entirely  different  meanings  of  d.  Therefore,  the  relative 


merits  of  integral  solutions,  with  and  without  the  present 
improvement,  oan  only  be  detei^ined  by  comparison  with  the 
corresponding  exact  solutions  *  In  this  regard^  it  is  clearly 
seen  from  fable  1  that  within  the  range  of  validity  of  the 
integral  solution  based  on  polynomial  profiles>  the  degrees 
of  accuracy  of  the  improved  solution  based  on  polynomials  far 
exceed  that  based  on  exponential  function.  All  these  observa* 
tions  are  further  substantiated  in  Pigs.  2-7 t  mci.,  which 
are  plots  of  f'^CT)^)  based  on  various  solutions  calculated 
for  different  a«values.  fhese  comparisons  of  f'''('n^)i 
which  are  directly  related  to  the  second  profile  derivatives, 
are  much  more  critical  than  that  of  f"(o)>  since  it  is  gen^ 
erally  kno%m  that  by  far  the  major  portion  of  inaccuracy  in 
the  basic  integral  solution  occurs  in  the  second  profile  de* 
rivatives  and  furthermore,  these  comparisons  are  made  over 
the  complete  range  of  Tlj^.  fhere  is,  however >  also  an  indi^ 
cation  that  the  accuracy  of  the  Improved  solution  decreased 
as  laminar  separation  is  approached* 

On  the  basis  of  the  above  coimDarisons  in  the  present 
problem,  which  does  cover  a  wide  range  of  free«stream  unsteadl 
ness,  it  is  perhaps  reasonable  to  suggest  that  in  a  general 
unsteady  problem  the  improved  integral  solution  based  on 
fourth»degree  velocity  pGiynomials  is  to  be  used  within  its 
entire  range  of  validity,  i.e.  I2^x^-12.  Outside  this  region 
for  X>12,  the  results  of  this  problem  suggests  the  use  of  the 
basic  Integral  solution  with  the  assumed  exponential  profile. 
However,  as  one  of  the  subsequent  examples  will  show,  it  is 
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not  always  tme  that  In  this  region  the  basic  integral  solu*- 
tion  is  more  accurate  than  the  present  solution.  Consequently^ 
it  is  necessary  in  this  region  to  evaluate  both  J*values  as 
to  determine  which  is  the  more  accurate.  CMe  obvious  wea^ess 
in  the  above  recommendation  is  the  use  of  composite  solution 
when  X  does  exceed  12.  However,  this  is  not  considered  to 
be  too  seriouSi  since  such  x*values  would  correspond  to  ex-^ 
tremely  large  accelerations  in  the  free  stream,  Purther 
justification  of  this  recommendation  can  be  obtained  in  con^ 
sidering  the  following  examples, 

ElCAMPLE  2,  The  Rayleigh  Problem  with  Step-Change  in  Freeg 
Stream  yelocity 

This  is  the  classical  problem  of  an  infinite  plate  moving 
in  a  stationary  incompressible  fluid  with  a  step-change  in 
its  velocity.  The  solution  is  well  icnown  and  is  given  by 
u*  s  erf(y/2/W)i  relative  to  the  coordinate  system  fixed  on 
the  plate,  where  erf  is  the  usual  error  function.  Solution 
to  the  integral  equation  in  this  case  is  elementary,  and  is 
shown  in  the  following: 
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The  integral  solution  Is  then  obtained  by  substituting  the 
above  equation  in  the  profile  equation  (10)  and  (12),  The 
corresponding  error  criterion  is  simply 
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where  ^  improved  integral 

solution  for  this  profelem  is  equally  elementary.  Punctions 
and  Pg  in  iquation  (29)  and  (30)  reduce  to 

^’lO  “r^  ^2Q  *  0 

respectively,  solution  to  iquation  (31),  satisfying  the 
boundary  conditions  in  {32),  is  simply  u*  =  erf(fl/$).  How'^ 
ever,  since  it  can  be  shown  readily  that  Tig  defined  as 
y/2yvt  is  related  to  J1  by  t)  =  Tig#,  the  solution  now  becomes 

u*  ^  erf  (Tig)  (4§) 

which  is  obviously  identically  to  the  exact  solution.  Thus, 
it  is  seen  that  for  this  problem,  the  improved  procedure  leads 
to  exact  solution,  regardless  of  the  type  of  profiles  used  in 
the  integral  solution.  The  calculated  J^values#  ns  well  as 
the  Surface  profile  derivatives,  are  shown  in  Table  2. 

EXAMPI^  3.  The^  Rayleigh  Problem  with  Constant  Acceleration 
This  problem  is  similar  to  the  previous  one  except  that 
the  free*stream  velocity  undergoes  constant  acceleration  from 
rest.  Exact  solution  to  this  problem  is  again  icnown  [6], 
and  is  given  by 


u*  =  1  +  ^  -  (1  +  [}  "  srf(Tig 

V  TT  ~ 
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.g/J  (^6) 

The  solution  to  the  integral  equation  (16)  has  been  given  by 


and  may  be  written  as 
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>  polj^&inlal 
Integral  Solution  | 

'  exponential 

Improved  Solution 
Exaet  Solution 


du*  J 

^dTlg  w  t* 


1.095^  3.4249 

1.2247  4.6340x10“^ 

1*1284  0 

1*1284  0 


vt 


B  H(K^) 


where  is  a  constat  satisfying  Iquation  (18).  Since  the 
GOrrespondlng  value  of  X  exceeds  +12,  only  the  result 
based  on  exponential  function  cah  be  considered*  Thus, 
Equation  (l8)  leads  to  a  devalue  of  "1/3 »  the  error  criterion 
for  the  Integral  solution  Is  now 
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Now  for  the  Improved  solution.  Equations 


and  (30)  reduce 


Pjo'IcC  '  l)  («9) 

The  integrals  in  the  solution  to  Equation  (31)  have  been 
evaluated  numerically,  yielding  the  new  profile  u*  =  u*('n). 
Finally,  this  may  be  transformed  into  the  physical  plane  by 
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Slftce  based  on  the  newly  obtained  profile  is  identical 

to  P^Q,  the  error  quantity  now  sifripiy  becomes  (u^  *  u  )/3 

in  view  of  Equation  (37)^  and  the  corresponding  J** value  has 
also  been  determined  from  Equation  (48) *  The  comparison  of 
J“ values  and  surface  profile  derivatives  is  shown  in  fable  3, 
and  that  of  the  second  profile  derivatives,  in  Fig*  8^  it 
is  seen  that  the  present  solution  is  much  more  accurate  than 
the  basic  integral  solution,  a  result  which  is  contrary  to 
that  obtained  in  the  first  example.  It  is  therefore  clear 
that  in  this  region  of  application  for  x>12>  it  is  necessary 
to  compute  the  error  quantities  for  both  the  integral  solutions 
with  and  without  the  present  improvement.  Such  that  the  more 
accurate  result  can  be  determined. 


TABLE  3  -  RiSULfS  FOR  EXAMPLE  3 


Integral  Solution  (exponential) 

du* 

'sr’w 

d-Hg 

2.3094 

J 

0.4459x10''- 

Present  Solution 

2.2662 

0.7682x10*^ 

Exact  Solution 

2.2568 

0 

Constant  AcGelerafelon  from  Rest 


This  fepresents  one  of  the  problems  in  which  the  leading 
edge  plays  a  very  Important  role  in  the  soliitionSi  These 
problems  have  been  systematically  studied  in  the  literature 
|^1^«14]  i  In  particular,  this  very  problem  is  covered  in  the 
series  solutions  of  Cheng  ^13} ,  and  Cheng  and  llliott  [14] , 
which  are,  however,  only  valid  for  flat  plates  and  free**-stream 
velocity  starting  from  rest*  In  this  numerical  example, 
attention  is  to  be  placed  only  in  the  region  close  to  the 
leading  edge,  since  it  is  in  this  region  that  the  basic  in¬ 
tegral  equation  (6)  for  this  problem  reduces  to  Equation  (17) 
exactly,  which  is  now  to  be  solved  in  terms  of  the  boundary 
condition  K(0)  >■  0*  The  Solution  may  be  eiscpressed  in  terms 
Of  integrals,  as  follows: 


which  immediately  yields  either  X  »  x(P)  or  Q  =  G(P), 
depending  on  the  assumed  profiles,  which  then  leads  to  the 
profile  details  according  to  the  integral  solution.  For  the 
improved  solution,  functions  P^^  and  Pg  now  become 
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f\.  X* 

-(X  «  2PX')  + 

2  2 
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+  XX'  «» 


ax  J 


u 


»  X  ^(u*  ^  2P)X'  i  +  u*  ij 


respectively,  where  prime  deftotes  derivatives  with  respect  to 
Pi  Iqiuations  (31)  are  explicitly  written  for  the  case  where 
velocity  polynomials  are  considered*  However >  they  are  equally 
valid  for  exponential  fimctions,  provided  again  that  x  is 
replaced  by  (1  +  2Q).  when  Iquations  (31)  are  evaluated 
from  the  integral  solution  Uq*  and  the  following  Identities 
are  utilized: 


X' 


2F3  3KH 
1  ^  3Pf| 


2(1  -  2a)(i  +  30) 

0(8  +  30  -  24P) 

where  prime  again  refers  to  derivatives  with  respect  to  P, 
Equation  (31)  IS  then  solved  in  ternis  of  integrals  which  In 
turn  are  evaluated  niunerically .  This  calculation  has  been 
carried  out  for  both  profiles  for  a  range  of  P  In  the 
nelshborhood  of  P  ^  Q,  The  J« values  for  the  basic  integral 
solutions  can  be  readily  evaluated  from  Equation  (36).  Mow= 
ever,  calculation  of  that  of  the  eorresponding  improved  solu-' 
tlons  gives  rise  to  a  slight  complication  which  needs  special 
treatment.  This  Is  a  characteristic  to  all  general  problems 
with  nons’Sirailar  profiles*  as  already  pointed  out  by  Yang  [27}. 


In  order  to  evaluate  the  error  quantity  a  from  Equation 
(37)i  it  is  neGessary  first  to  determine  functions  and 
Pg  from  the  improved  profiles  *  However,  in  view  of  Equations 
(51),  the  distribution  of  3uV^P  at  all  values  of  T],  must 
be  known.  Even  though  it  could  possibly  be  evaluated  by 
numerical  differentiation  from  improved  profiles  at  various 
P- values,  such  an  evaluation  would  be  extremely  awkward,  and 
it  is  highly  desirable  to  be  able  to  determine  this  distrlbu* 
tion  at  any  particular  value  of  P.  This  c^  be  accomplished 
by  the  following  scheme*  Different latlng  Equation  (3I)  with 
respect  to  P  yields 

aa  Su*  a  a, go  s,,.  au* 

_(— )  +  p  — ~  -  5.^2  —  (sa) 

dT|2  dp  -  ^  dT|  dp  dp  dp  dT] 

which  is  subjected  to  the  boundary  conditions 
du*  du* 

T1  »  0  ^-0  ^^0  (53) 

dp  dp 

Now  Equation  (52)  may  be  solved  in  an  identical  way  as  that 
for  Equation  (31),  resulting  in  du*/dP.  Thus,  this  determina 
tion  is  frozen  at  any  p- value. 

The  calculated  results  for  this  problem  are  shown  in 
Pigs.  9  and  10,  together  with  the  series  solution  of  Cheng 
^13] ,  which  is  believed  to  be  accurate  In  this  problem,  in 
view  of  the  nature  of  his  solution.  More  specifically.  Pig. 

9  shows  the  variation  of  local  skin*frlction  coefficient 
Cf  -  T^/(l/2pUp~)  with  p  and  Pig.  10,  the  variation  of  the 
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J= value Si  It  Is  noted  that  all  curves  in  Pig*  10  go  through 
J  =  0  at  P  »  Oi  slnee  there  is  no  motion  then  and  accord* 
ingly,  6  =  0,  yielding  j  =  0  in  view  of  Equation  (36) i 
fhese  results  again  show  the  good  correlation  hetween  the 
calculated  values  and  the  relative  accuracies  of  the  integral 
solutions  with  and  without  the  present  improvement,  this 
represents  further  evidence  of  the  validity  of  the  use  of  the 
error  criterion.  Moreover >  it  is  also  seen  that  the  improved 
integral  solution  based  on  polynomials  is  much  more  accurate 
than  that  based  on  exponential  functions.  Since  in  this 
problem  values  of  x  are  all  under  the  limit  of  12 >  this  re* 
suit  agrees  well  with  that  of  the  first  example. 

EXiUMPIiE  |.  f ransienfe  Boinidary*I«ayer  Development  on.  a  Settl* 
Infinite  Plate 

This  problem  deals  with  the  effect  of  step*wlse  change  in 
free*stream  velocity  on  the  development  of  unsteady  boundary 
layer  on  a  eeml-infinite  plate  which  is  originally  in  steady 
motion.  The  step*wise  change  may  represent  either  accelera* 
tion  or  deceleration.  The  integral  solution  to  this  problem^ 
based  on  the  method  of  characteristics «  has  already  been 
presented  previously,  and  is  shown  in  Equation  (23).  Only 
the  transient  solution  is  considered  here,  since  the  steady'^ 
state  solution  has  already  been  treated  in  the  literature 
^26l.  For  this  problem,  functions  and  ?2  4re  given  by 
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Detailed  ftitmefical  GalGulations  of  the  integral  solutions 
with  and  without  the  present  Improvement  >  hased  on  both  as^ 
sumed  profiles >  have  been  Garried  out  for  two^vaiues  of 

namely,  2*0  and  0,5,  and  several  eomblnations  of  x* 
and  t*»  In  addition>  all  pertinent  values  of  J  have  also 
been  determined*  since  n©  exaGt  solution  is  available  to 
this  problem  in  the  literature,  no  Gomparlson  ean  be  made 
and  only  some  representative  results  are  presented*  fable 
4  shows  the  GalGulated  results  for  several  speGifle  oombina^ 
tions  of  X*  and  t*,  and  some  typioal  velocity  profiles  are 
described  in  Pig,  11  to  14,  incl*  fhese  results  again  clearly 
indicate  the  validity  of  the  improved  solution.  Pinally>  it 
is  to  be  noted  that  in  a  normal  application  to  the  present 
problem,  it  is  not  necessary  to  obtain  the  integral  solutions 
based  on  exponential  functions >  since  here  x  is  identically 
zero,  which  is  well  within  the  range  of  validity  of  the  in* 
tegrai  solutions  based  on  polynomials. 
example  6 .  tinst^eady  Plow  over  a  Circular  Cylinder 

This  last  exs^le  treats  a  general  unsteady  boundary’* 
layer  problem  of  flow  over  a  circular  cylinder*  The  free** 
stresum  velocity  distribution  is  written  as 

=  d(x*)s(t*)  (55) 


0(x*)  =  3.63i4x*  -  2.l709x*^  **  l.5l44x*^ 

g(t*)  *  1  *  t*  +  t*^  t*^Q 
*  I  t*<0 


where 
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where  0(x*)  is  the  Steady  free^stream  velocity  distributioh 
over  a  circular  cylinder  as  given  by  Heimenz  [32] .  Even  though 
this  variation  is  now  known  not  to  be  too  accurate,  it,  never* 
theless,  serves  the  present  purpose.  The  chosen  tiine  varia* 
tion  is  one  involving  initially  steady  motion,  then  decelera* 
tion,  and  finally  followed  by  an  acceleration «  The  general 
calculation  procedure  described  previously  for  a  general  un* 
steady  problem  to  obtain  solution  to  the  integral  equation  is 
here  used.  All  integrations  have  been  carried  out  by  the 
usual  fourth-order  Runge-Kutta  routine.  The  result  of  this 
solution  based  on  poljmomial  profiles  only  is  shown  in  Pig.  13 » 
For  this  problem  Equations  (29)  and  (30)  are  used  directly, 
Since  no  simplification  is  possible.  When  the  integral  soiu- 
tion  Uq  (XiTl)  is  substituted  into  Equations  (33)>  improved 
profile  at  any  combination  of  x*  and  t*  is  easily  obtained. 
To  calculate  the  error  quantity  of  the  improved  solution,  it 
is  now  necessary  to  first  evaluate  duV^x  and  duV^t* 
from  the  improved  profile.  This  is  again  done  by  first  dif* 
ferentiating  Equation  (31)  with  respect  to  \  and  t*,  and 
then  following  the  same  scheme  as  that  indicated  in  iquation 
(32).  Some  results  are  shown  in  TabDe  9  and  Fig.  16.  Here 
again  the  improved  accuracy  in  the  present  solution  may  be 
noted. 


TABLE  5  * 


t* 

X* 

G.  10000 

0,20000 

0 . 30000 

0,10000 

0.20000 

0.10000 

0.10000 

0.40000 

o.ioooo 

0,40000 

0.30000 

0.30000 

0,40000 

0. 20000 

0.40000 

0.40000 

0.30000 

0.30000 

0.50000 

0.40000 
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CONCLUDING  REMARKS 

In  this  report  an  improved  integral  prOGedure>  based  on 
two  types  of  assumed  profiles,  is  proposed  to  calculate  the 
behavior  of  unsteady  laminar  boundary  layers  over  arbitrary 
cylinders  with  arbitrarily  prescribed  unsteadiness  in  the 
free  stream*  Also  introduced  is  a  simple  error  criterion 
by  which  the  validity  of  the  Improved  solution  can  be  readily 
determined.  In  view  of  results  from  the  numerical  examples, 
this  error  criterion  correlates  extremely  well  with  the  in^ 
accuracy  of  the  approximate  solutions  on  the  basis  of  comparing 
results  with  that  from  the  exact  solutions.  Furthermore,  the 
critical  comparisons  in  the  second  profile  derivatives  from 
the  Integral  solutions  with  and  without  the  present  improve* 
ment  technique  and  from  the  exact  solutions  in  the  first  four 
nurnerlcal  examples  have  definitely  indicated  the  high  degree 
of  accuracy  attainable  in  the  present  solution,  as  mentioned 
previously,  such  accuracy  is  necessary  for  hydrodynamic  stab¬ 
ility  considerations.  Consequently,  the  following  recoiranenda- 
tion  has  evolved  in  the  present  study.  For  any  general  unsteady- 
flow  problem,  the  improved  integral  solution  based  on  the 
fourth-degree  polynomial  as  the  assigned  profiles  in  the  basic 
integral  solution  should  be  used  within  the  entire  range  of 
validity  of  this  basic  integral  solution.  High  degree  of 
accuracy  c^  be  expected,  except  possibly  in  the  immediate 
neighborhood  of  the  point  of  laminar  separation.  Outside  of 
this  range,  both  integral  solutions  based  on  asstmied  expon* 
ential  functions  should  be  utilized.  The  final  solution  here, 
being  either  the  basic  integral  solution  or  the  improved  solution, 
is  evidently  the  one  with  lower  magnitude  of  the  error  quantity  J. 
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